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Abstract. 

Relativistic phenomena (orbital precession, Shapiro delay, and/or 
orbital decay) have been measured in Arecibo timing observations of four 
pulsar-white dwarf binaries, leading to constraints on the neutron star 
masses. We have detected the decay of the PSR J0751+1807 orbit due to 
gravitational radiation emission, the first such measurement in a binary 
with a low mass ratio {m2/mi ^1). Timing data constrains the mass of 
this pulsar to bet between 1.6 and 2.8 Mq. Masses of the other pulsars 
are in marginal agreement with the canonical pulsar mass of 1.35 Mq, 
but higher values are preferred. 



1. Introduction 

The conventional wisdom is that all pulsars have mass close to the Chan- 
drasekhar value, ~ 1.35 Mq. This belief has arisen from high precision measure- 
ments of neutron star-neutron star binary systems; in all cases, the pulsar and 
the secondary neutron star have masses within a narrow range, 1.33 to 1.45 Mq 
(Stairs et al. 2002, Weisberg & Taylor 2003, Deich & Kulkarni 1996). Perhaps 
all neutron stars are born within this narrow mass range, or perhaps the particu- 
lar evolution of double neutron star leads to this result (Brown 1995). Even if all 
neutron stars are born near the Chandrasekhar mass, they may achieve higher 
masses through accretion of matter from orbital companions. Millisecond pul- 
sars in orbits with white dwarfs, in particular, are known to have accreted mass 
during their low-mass X-ray binary phases, and hence they are prime candidates 
to be relatively massive pulsars. There has been some observational evidence for 
heavy pulsars in these systems (Kaspi, Taylor, & Ryba 1994; van Straten et al. 
2001), although a statistical study argued that they are no heavier than pulsars 
in double neutron star binaries (Thorsett & Chakrabarty 1999). 

In this status report, we summarize observations of four millisecond pulsars 
in pulsar-white dwarf binaries made with the Arecibo telescope. This is an 
update of Nice, Splaver, & Stairs (2003), incorporating significant new timing 
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Table 1. Mass Measurements from Timing Analysis 
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data from PSR J0751+1807 and an improved timing model (including orbital- 
annual parallax) for PSR J1713+0747. 

2. Observations &; Data Reduction 

Observations were made at Arecibo over several years at radio frequencies 430 
and 1410 MHz, primarily using the Princeton Mark IV data acquisition system 
(Stairs et al. 2000). Details will be given elsewhere. These data were augmented 
by observations at Green Bank (140 Foot telescope), Jodrell Bank, and Effels- 
bcrg, as well as data collected at Arecibo using the Arecibo Berkeley Pulsar 
Processor and the Princeton Mark III system. 

Pulse arrival times were derived using conventional methods and were fit 
to timing models using TEMPO^. For each pulsar binary, at least one post- 
Keplerian orbital parameter is detectable, leading to constraints on pulsar mass, 
nil, secondary star mass, m2, and orbital inclination, i. (See Taylor, 1992, for the 
relations between post-Keplerian parameters, masses, and orbital inclination.) 
To derive confidence intervals on mi, m2, and i, we calculated the goodness-of- 
fit of timing solutions over a grid of trial values of cosz and m2- Each pair of 
trial values corresponds to a particular pulsar mass, mi, according to the mass 
function, fi = (m2sini)^(mi +m2)~^ = {2Tr / Pf,)'^ /Tq, where Tq = 4.925 /xs 
is the solar mass in time units, Pj, is the orbital period, x is the semi-major axis 
projected into the line-of-sight in time units, and mi and m2 are in solar mass 
units. For each trial, relativistic post-Keplerian parameters appropriate to cosi 
and m2 were calculated and held constant, while all other parameters were fit. 

The contours plotted in Figure 1 show the combinations of cos i and mi that 
gave acceptable fits for each source. We expect cos i to be uniformly distributed 
for randomly oriented binary systems. (Other than this, the intrinsic values of 
i are of limited interest.) Table 1 summarizes the pulsar parameters and 95% 
confidence limits on pulsar masses from the timing data. 

3. Individual Sources 

PSR J0621+1002 . This is an "intermediate mass" binary millisecond pulsar, 
presumably with a CO white dwarf secondary. Constraints on the masses come 



^ http: / /pulsar. princeton.edu/tempo 
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Figure 1. Constraints on pulsar masses and orbital inclinations. In- 
ner and outer contours are 68% and 95% confidence regions, respec- 
tively, based on timing data alone. Shaded areas indicate theoretical 
limits within these contours, from the assumption that the companion 
is a white dwarf (J0621+1002) and from predictions of secondary star 
masses from the Pb-m2 relation (B1855-I-09 and J1713+0747). (There 
are no theory constraints on J0751+1807.) 

from the relativistic advance of periastron — we measure the precession rate to 
be cj = 0.0116 ± 0.0008° yr"^ {la), yielding mi+m2 = 2.8 ± 0.3 Mq— and from 
the lack of detection of Shapiro delay, which excludes large inclinations. While 
the canonical pulsar mass value of 1.35 Mq falls within the allowed region (Fig. 
1), higher values are preferred, especially when we include the constraint that 
the white dwarf secondary have mass less than 1.4 Mq (shading in Fig. 1). See 
Splaver et al. (2002) for more details. 

PSR J0751 + 1807. This pulsar has a very short orbital period. By com- 
bining pre- upgrade (1993-1994) and post-upgrade (1999-2003) Arccibo obser- 
vations of this source, along with data from Jodrell Bank and Effelsberg, we 
are able to detect the orbital decay due to emission of gravitational radiation. 
The rate of change of orbital period is Pf) = (—6.2 it 1.1) x 10~^^. This is the 
first detection of the decay of a pulsar-white dwarf orbit with a low mass ratio 
(m2/mi <^ 1).^ Shapiro delay is marginally detected, constraining i to interme- 



"^Notc also the spectacular results from PSR J1141— 6545, a young pulsar in an eccentric orbit 
with a massive white dwarf (Bailes et al. 2003.) 
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diate values. If the Shapiro delay detection can be confirmed and refined it will 
provide independent constraints on the system masses, allowing for novel tests 
of relativistic gravity (Arzoumanian 2003). In any case, the allowed values of 
nil are significantly greater than 1.35 Mq (Fig. 1). 

PSRs B 1855 +09 & J 171 3 +0747 . These pulsars are in wide orbits with hght 
white dwarf companions. In both cases, the Shapiro delay is detectable, con- 
straining i and mi as shown in Figure 1. Binary evolution theory posits a narrow 
range of companion masses, m2, for a given orbital period, Pb, in these systems 
(Podsiadlowski, Rappaport, & Pfahl 2002; Tauris &: Savonije 1999; Rappaport 
et al. 1995). The resulting constraints on mi and i are indicated in Figure 1. 
Once again, values of mi somewhat higher than 1.35 Mq are preferred.^ 

4. Conclusion 

The sources individually provide some evidence for pulsars substantially more 
massive than the Chandrasekhar value, and the collection of constraints shown 
in Figure 1 provides a compelling case for heavy neutron stars. 
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^Parkes pulsar J0437-4715 gives a similar result, mi = 1.58±0.18 Mq (van Straten et al. 2001). 



